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Abstract 

We consider the Ising model on an M x rectangular lattice with an asym- 
metric self-dual boundary condition, and derive a closed-form expression for 
its partition function. We show that zeroes of the partition function are given 
by the roots of a polynomial equation of degree 2M — 1, which trace out cer- 
tain loci in the complex temperature plane. Particularly, it is shown that (a) 
real solutions of the polynomial equations always lead to zeroes on the unit 
circle and a segment of the negative real axis, and (b) all temperature zeroes 
lie on two circles in the limit of M — > oo for any A^. Closed- form expressions 
of the loci as well as the density of zero distributions in the limit of — > oo 
are derived for M = 1 and 2. In addition, we explain the reason of, and es- 
tablish the criterion for, partition function zeroes of any self-dual spin model 
to reside precisely on the unit circle. This elucidates a recent finding in the 
case of the self-dual Potts model. 
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Typeset using REVT^ 



I. INTRODUCTION 



One of the most remarkable results in statistical mechanics is the circle theorem of Yang 
and Lee who in 1952 established that zeroes of the partition function of the ferromagnetic 
Ising model always lie on the unit circle in the complex e~^^ plane, where L = H/kT is the 
reduced magnetic field variable. In 1964 Fisher p[ pointed out that it is also meaningful to 
consider zeroes of the Ising partition function in the complex temperature plane, in which the 
zeroes reside on a boundary at which the free energy becomes non-analytic. Particularly, 
he showed that temperature zeroes of the square lattice Ising model lie on two circles in 
the thermodynamic limit. Since then the consideration of the partition function zeroes has 
become a power tool in analyzing lattice models. The temperature zeroes of the Ising model 
1^,^ as well as the Potts model have been investigated. However, studies to this date have 
been confined mostly to results in the thermodynamic limit and/or numerical analyses 
I^J^. Very little progress has been made for finite lattices, except in one dimension [Q, in 
which the locations of zeroes are determined in a closed-form and analyzed algebraically. 

In this paper we study the partition function zeroes of the Ising model for a finite rect- 
angular lattice with an asymmetric self-dual boundary condition. This study is motivated 
by a recently discovered numerical evidence in the case of self-dual Potts models |^ that 
many partition function zeroes lie precisely on the unit circle. Here, we explain why this 
is the case, and establish more generally the criterion for zeroes to occur on the unit circle 
for self-dual lattice models. Particularly, for the Ising model considered here, we show that 
all zeroes lie on two circles in the limit of M — > cxd and any A^. We also obtain closed-form 
expressions for the loci of zeroes for M = 1, 2, as well as the density of zero distribution in 
the limit of iV — >• oo. 



II. A SELF-DUAL RECTANGULAR LATTICE 

Consider a rectangular Ising lattice C of M rows and columns with periodic boundary 
conditions in the horizontal direction. Let the nearest-neighbor interaction be —kTK. Fur- 
ther, introduce a boundary condition that all sites on one of the two horizontal boundaries, 
say, the lowermost one, interact with an additional site with the same interaction —kTK. 
Thus, there are altogether MN + 1 sites and 2MN edges. Note that the boundary condition 
is "asymmetric" in the horizontal and vertical directions. 

Topologically, £ assumes the form of a "wheel" consisting of N spokes and M concentrate 
circles. An example of C is shown in fig. 1(a) where the circumference of circle corresponds 
to the "horizontal" direction. We point out that the lattice C is self-dual, which is an 
important premise of our consideration. 

The high-temperature expansion of the partition function assumes the form 

Zm,n = 2^^^+i(coshir)2^^^Ga^) (1) 

where z = tanh K and 

Gc{z) = l+ E (2) 
c.p.e c 
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Here, the summation in (0) is over all closed polygonal configurations that can be drawn on 
£, and b is the number of edges contained in each polygonal configuration. 

To facilitate our considerations, we expand the center point of the wheel in Fig. 1(a) 
into a circle as shown in Fig. 1(b), and associate edge weights 1 to the M newly added 
edges. This transforms the lattice C into C*. A moment's reflection [j^ shows that we have 
the relation 



(3) 



It follows that the partition function (P is computed if we can evaluate G*^{z). 

The lattice C* is an (M + 1) x rectangular lattice with a periodic boundary condition 
in the horizontal direction. All edges of C* carry the weight z except those on the innermost 
circle which carry the weight 1. Following Kasteleyn [§] we can express Gc*{z) as a Pfaffian 
which is, in turn, the square root of a determinant. The procedure, which is standard, has 
been described in details by McCoy and Wu P,|lO|, and leads to 



r N 



GrAz) 



ndet|5M(0)| 



n=l 



1/2 



where Bm{(P) is a 4(M + 1) x 4(M + 1) matrix given by 



/ B{z) B+{z) 
B.{z) B{z) B+{z) 



Bm{<P) 



V 



Here, the product in (^) is taken over the A^ values 



= (2n - l)7r/A^, 



(4) 



B_{z) B{z) B^{z] 
B-{z) B{1) ) 



(5) 



n 



1,2,- ■ -AT, 



(6) 



and B{z) and Bjjyz) are 4x4 matrices given by ( [A5| ) in the Appendix. For completeness, 
an outline of the derivation of (|^) is given in the Appendix. Note that 7^ 0, vr. The 
substitution of (^ into (|^) and (|I]) now expresses the partition function Zm,n in terms of a 
4(M + 1) X 4(M + 1) determinant. 



III. EVALUATION OF THE PARTITION FUNCTION 

In this section we evaluate the determinant det \Bm{<^^\ in (^) using an approach some- 
what different from that of P,p^, and present the solution in a form more suitable for our 



purposes. 

For brevity we write Bm = det \Bm{4')\- The determinant Bm is bilinear in the two 
non- vanishing elements z in B^{z) in the first row and —z in B„{z) in the first column of 
(§). Explicitly, the bilinear expansion is 

Bm = det \B{z)\ ■ Bm-i + z[BM]ifi - z[BM]fi;i + z^[Bm]i,6;1,6 (7) 
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where [BM]i;k is the determinant of the matrix BM{<f>) with the ith row and kth column 
removed, and [BM]{ij-ke) is similarly defined. We have the following lemma whose proof is 
elementary and will not be given: 

Lemma: For any any m x n matrix Amn and (A^ — m) x (A^ — n) matrix CN-m,N-n-, we 
have the identity 



det 



A n 

CN~m,N~n 



Sm,n ■ det \Amn\ ' det \CN-m,N-n\, (8) 



where 6 is the Kronecker delta function. 
It follows from the lemma that we have 

[BMh-fi = [Bm]6;1 = 0, (9) 
[5Af]l,6;l,6 = det • [5m-i]2;2. (lO) 

Writing Dm = [Bm]2,2 and evaluating det |-B(2;)| and [-Ba/]i,6;1,6 explicitly, (|^) becomes 

Bm = aBM-i + bDM-u (H) 

with 

a = 1 + z'^ - 2zcos(j), b = 2iz^sm(j). (12) 
In a similar fashion we expand Dm and obtain 

Dm = cBm-1 + dDM-i, (13) 

with 

c=-2izsm(j), d = z^{l + z'^ + 2zcos(j)). (14) 
Both (|l^) and (|T3D hold for any M > 1 with the boundary condition 



5o = det|5(l)| =2-2cos(/), Do = det |5o(l)|2,2 = -2isin(/). (15) 
Introducing generating functions 

oo oo 

B{t) = Y: Smt"", D{t) = Y: Omt"" (16) 

Af=0 A/=0 

and multiplying both sides of ( |Tl| ) and (|13[) by r^, followed by summing over from M = 1 
to oo, one obtains 

B{t) -Bo = T[aB{T) + bD{r)] 

D{t) -Do = t[cB{t) + dD{T)]. (17) 

This yields 
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Bir) 



Bo{l-dT) + Dobr 
1 — (a + d)T + {ad — hcjr'^ 



(18) 



It is then a simple matter to obtain an explicit expression for Em by expanding the right- 
hand side of (0) and comparing with (p!6D, after first carrying out a partial fraction. After 
some algebra, this leads to the expression 



2\iM 



Bm = {2~2cos(j))[z{l-z')] 



sinh[(M + l)t(0)] - s-^ sinh[Mt(0)] 



sinht(0) 



where 



2 cosht(0) = s + s ^ — y 

s = s{K) = e^^ cothK, |/ = 2(l + cos< 



(19) 



(20) 
(21) 



Substituting (|T9D into (0) and (^, we then obtain from ([l|) the following closed- form expres- 
sion for the partition function: 



n=l 



sinht(0) 



(22) 



Here, use has been made of the identity |]TT 



N 

n 

ra=l 



2 - 2 COS 



/2n - 1 



N 



-TT 



4. 



(23) 



We remark that we find also 



D 



M 



2mA/ 



sm I 



sinht(0) 



sinh[(M + l)t(0)] + s{-K) sinh[Mt(0)] 



(24) 



IV. ZEROES OF THE PARTITION FUNCTION 

/^From the closed-form expression (^21) we find the partition function to vanish at 

^ sinh[Mt(0)] 

sinh[(M + l)t(0)]' ^ ^ 

a result which holds for all M > 1 and equal to the N values given by (P). Note 
that, since only cos0 appears in the expression, distinct zeroes are obtained by using only 
n = 1, 2, ■ ■ ■ , [(A^ + l)/2], where [(A^ + l)/2] is the integer part of (A^ + l)/2, equal to A^/2 
for A^ = even and (A^ + l)/2 for n = odd. It should also be noted that the numerator and 
denominator of (pSD contain a common factor sinht(0) which cancels, resulting in polyno- 
mials of 2cosht(0), or s + — y, of degrees M — 1 and M, respectively. Equating both 
sides of (p5D after multiplying by the denominator, one finds that the two terms of the order 
of s~^^ cancel, so that ( [25| ) assumes the form 
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P2M-iis) = 0, n = l,2,---,^^— 

where P2m-i{s) is a polynomial of degree 2M — 1 in s. Explicitly, we find 



(26) 



2cosht(0) - s"^ = 0, M = l 

4cosh2t(0) -2s~^cosht(0) - 1 = 0, M = 2 

8cosh3t(0) -4s"^cosh2t(0) -4cosht(0) + s"^ = 0, M = 3, (27) 



which yield, after using (^ 



s - y = 0, M = 1 (28) 

- 2ys^ + y^s-y = 0, M = 2 (29) 

s^-3ys'^ + 3y^s^-{2y + y^)s^ + 2y^s-y = 0, M = 3. (30) 

Note that no closed- form expression exists for the roots for M > 3. 

Guided by results of a recent study 0, we now introduce the variable 



X = (e^^ - 1)/V2 (31) 

so that from (^Tj) we have 

s = V2{x + x-^) + 3. (32) 

This says that if x is a root of (|26|) , then is also a root. Furthermore, since |a; + a;^^| < 2 
if, and only if, x is on the unit circle |x| = 1, all real solutions of (p6|) in the regime 



3 - 2V2 < s < 3 + 2V2 (33) 

lead to solutions on |x| = 1. In fact, more generally, we have the following theorem: 

Theorem 1: Let Pl{x) be a polynomial of degree L in x satisfying the reciprocal relation 

Pl{x)=x'^Pl{x-'). (34) 

Then we have (i) Pl{x) can be rewritten as a polynomial in a new variable w with 

Pl{x) = x^^'^Ql/2{w), w = X + x^^, L = even 

= x^/^Ql{w), w = x^/^ + x~^/\ L = odd, (35) 

where Ql{w) is a polynomial of degree L in w. 

(ii) Roots of Pl{x) = are on the unit circle |a;| = 1 if, and only if, the corresponding 
roots of the polynomial Q{w) = are real and in the regime —2 < w <2. 
Proof: 

Property (i) follows from the elementary facts that (a) the reciprocal relation (|3^ implies 
X and x~^ appear symmetrically in the Q polynomials in the combination of X„ = x^ + x~^, 
where n is an integer for L even and a half-odd integer for L odd, and (b) X„ is a polynomial 
of w. Property (ii) follows immediately. Q.E.D. 
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A consequence of Theorem 1 is that zeroes of the partition function of self-dual spin 
models, which by definition satisfy the reciprocal relation ( |34D in an appropriate dual variable 
X, reside precisely on the unit circle |x| = 1 for all real roots of the Q polynomial in \w\ < 2. 
This explains the numerical findings of that a large number of zeroes of self-dual Potts 
partition functions are located on the unit circle. 

For the self-dual Ising model on C, we have further: 

Theorem 2: (i) All real roots of the [(A^ + l)/2] polynomial equations P2M-i(s) = are 
bounded by 

- 1 < Smin < S < Smax < 3 + 2v^, (36) 

where the bounds Smin and Smax depend on M. 

(ii) Zeroes in the x-plane corresponding to real s lie either on the negative real axis 
-(v^+1) <x_(s^i„) <x<x+(s^in) < -(v^-1) <0, s^in < s < 3 - 2^2, (37) 

or on the horseshoe segment 9q to 27r — 9q on the unit circle 

|x| = l, 3-2y2<s<w. (38) 

Here, Xj^{s)x_{s) = 1 with 

x±{s) = (s - 3 ± Vs^ - 6s + l)/2V2, 

^^o = cos-M(w-3)/2x/2]. (39) 

Note that the horseshoe segment becomes a circle (^^o = 0) if s^ax = 3 + 2-\/2. 

(iii) Zeroes in the x-plane corresponding to |s| = 1 lie on the circle \x + \/2\ = 1. 
Proof: 

(i) We first show that no root exists for s > sq = 3 + 2\^. Assume the contrary, namely, 
a solution s > Sq exists. Then by (|^) we have 

COsht(0) = -{s + S^^) — (1 + COS0) 

> ^(-50 + Sq^) - (1 + COS0) 

= 2 + cos 

> 1. (40) 

It follow that t{(f)) is real and ( ^5]) gives s < 1 which contradicts the assumption of s > sq = 
3 + 2 1/2. Hence no real root exists for s > sq. 
Likewise, for s < solution of fEBI) we deduce 



cosht(0) < -1 (41) 
for which t(0) must assume the form 

t(0) = r + i{2n + l)7r, (42) 
where r is the real part and n is an integer. Then (|25|) leads to 
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— sinh Mr 
sinh(M + l)r' 



(43) 



showing that negative solutions of (|26|) always lie in the interval (—1,0). This completes the 
proof of (i). 

(ii) For —1 < Smin < s < 3 — 2^/2, one solves x = x±{s) from (^). This leads to 
where we have used x±(— 1) = — (-\/2 =f 1). For 3 — 2^/2 < s < Smax, one writes x = e*^ 
solves for 6 from (|32D . This leads to the result (|38|) with 6q given by (0). 

(iii) It is straightforward to show that |s| = 1 by substituting x = —\/2 + e^^ into 
for any real (3. Q.E.D. 



and 



V. LOCI OF ZEROES 



The loci of zeroes can be explicitly determined for M 
Consider first the case of M = 1. 
For M = 1 ^ 



gives 



s = 2(1 + cos< 



1, 2 and in the limit of M — > oo. 



(44) 



so that Smin = and s^ax = 4. It follows that the zeroes are either on the line segment 
X = [—V2, —l/\/2] or on the horseshoe segment of = 1 with 6*0 = cos~"'^(l/2v^) = 57.3°. 
A plot of the loci obtained for = 4, 000 is shown in Fig. 2. In the limit of large A^, the 
zeroes are distributed uniformly in with a density p{(f)) = 1/2tt. Combining (^41), (^) and 
X = e^^, we have 



cos(f) = v^cos^ — 1/2, 
and find that the zeroes distribute on the horseshoe segment with a density 



1 sin 6' 



2tt dO v^TT sin i 



(45) 



(46) 



which diverges at 6q as {6 — 6q) 



For M = 2 the cubic equation (^) yields three roots 



22/ 



+ c^"VQ + v^ + ^'"VQ- n = 0,l,2. 



(47) 



where 



Q 



27 



A 



_ V 

4 V 27 



(48) 



Now So is always real and assumes the minimum value s = at y = 0. In addition, si^2 are 
real for A < or 3-\/3/2 < y < A. It follows that we have Smin = 0. A more detailed analysis 
shows that we have Smax = 8(1 + cosq;)/3, where a = | cos^^(— 5/32). These lead to zeroes 
in the x-plane on the negative real segment [— -\/2, — l/\/2] and a horseshoe segment of the 
unit circle from 6q to 27r — with 6q = cos^^[(8 cosa — l)/6-\/2] = 47.7°. A plot of the 
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zeroes for N = 2, 000 is shown in Fig. 3. The density of zeroes in the hmit of ^ oo, 
which can be found accordingly, is given as the superposition of the densities of the three 
branches. However, since only one of the overlapping densities extends to the density 
diverges at Oq as {9 — Oq)~^/'^. 

For M > 3 the polynomial equation (|2^) can be solved only numerically. We have carried 
out its numerical solution for 3 > M > 8 and determined the corresponding values of ^o- 
The result of M = 3, = 400 is shown in Fig. 4. We also show in Fig. 5 the trend of 
as M oo. 

In the limit of M — > oo, we note that for M large we must have le*'-'^-'! = 1. This is 
so since, otherwise, without the loss of generality we can assume le**^"^^! > 1 and ( pSf ) yields 
s ~ e"*^*^-* which contradicts (^) (except at = vr). Thus we write = ia{(j)) where a(0) 
is real. Then, (^) becomes 

s + = 2cosa + 2(1 + COS0) = real (49) 

leading to the bound —2 < s + s^^ < 6. For —2 < s + < 2 we must have |s| = 1 which 
implies, by Theorem 2(iii), that all zeroes are on the circle |a: + \/2| = 1- For 2 < s + s~^ < 6, 
the bound is equivalent to the bound (|33D on s and, by Theorem l(ii), all zeroes are on the 
unit circle |x| = 1. We therefore conclude that in the limit of M — * oo all zeroes are on 
the two circles |x| = 1 and \x + \/2\ = 1. Indeed, this coincides with the finding of Fisher 
who showed that, in the limit of both M, iV — >■ cxo, the partition function zeroes are on 
the two equivalent circles | tanhK ± 1| = |]l3l- However, our result is more general and 
holds for M —>■ oo and any A^. 

Finally, we remark that our results confirm the conjecture put forth in that all zeroes 
in the Re (x) > half plane are on the unit circle |x| = 1. But a generally proof of this 
conjecture is still lacking at this point. 

VI. DISCUSSIONS 

We have obtained the closed-form expression of the partition function of an M x 
rectangular Ising lattice with an asymmetric boundary condition, and evaluated its zeroes. 
It was found that the zeroes are given by the roots of the polynomial equation (^). Partic- 
ularly, we established that real roots of (|26D lead to zeroes either on the unit circle |x| = 1 
or within the segment [—(a/2 -|- 1), —{\/2 — 1)] of the negative real axis in the complex x 
plane. We also established that all zeroes lie on two unit circles centered respectively at 
X = 0, —V2 in the M ^ oo limit for any A^, where x = {e^^ — l)/y/2. Numerical results for 
finite M are also presented. 
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APPENDIX. EVALUATION OF A PFAFFIAN 



For completeness we include in this appendix an outline of the steps leading to expressing 
the graph generating function Gc*{z) as a Pfaffian. First, following standard procedures |§] 
we see that the graph generating function Gc*{z) for C* is equal to the dimer generating 
function for a "bathroom-tile" lattice constructed from C* as shown in Fig. 6. This leads 
to the relation 



Gc* = Pf A = ^det \Al 
where ViA is the Pfaffian of a 4(M + 1)A^ x 4(M + 1)A^ antisymmetric matrix 



Hn — C 



T 
M 



H 



A = Am O /iv + Cm 
Here, In is the N x N identity matrix, is an N x N matrix given by 



(Al) 



(A2) 



H 



N 



/ 
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• ^\ 
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. 
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-1 
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(A3) 



Am and Gm are4(M+l) x4(M+l) matrices given by 
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M 
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/ 5(0) 





G{z) 
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M 



B{0) 



G{z) 


Biz) 
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(A4) 
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o) 









(A5) 



Here Hj^ and Gj^ denote, respectively, the transpose of H^^ and Gm- 

Now, if AT and Hj^ can be diagonalized simultaneously with respective eigenvalues e^'^ and 
e~*"^, where = 7r(2n — 1)/A^, n = 1,2, ■ ■ ■ , N. We can therefore replace if at (resp. if^) by a 
diagonal matrix with diagonal elements e*"^ ( resp. e"*"^) without affecting the determinant. 
The substitution of (|A2D into (|AT|) now leads to M). 
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Figure Captions 



Fig. 


1. 


(a) A 3 X 8 self-dual lattice L. (b) The 4x8 lattice £* derived from L 


Fig. 


2. 


Loci of zeroes in the x plane for M = 1, N — 4:, 000. 


Fig. 


3. 


Loci of zeroes in the x plane for M = 2, N = 2, 000. 


Fig. 


4. 


Loci of zeroes in the x plane for M — 3, N — 400. 


Fig. 


5. 


The dependence of on M. 


Fig. 


6. 


The bathroom-tile lattice for a 3 x 4 lattice £*. 
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